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ABSTRACT 



This paper investigates which constraints can be placed on the fraction of Compton-thick AGN in the Universe from 
modelling the spectrum of the diffuse X-ray background (XRB). We present a model for the synthesis of the XRB 
that uses as input a library of AGN X-ray spectra generated by Monte Carlo simulations. This is essential to account 
for the Compton scattering of X-ray photons in a dense medium and the impact of that process on the spectra of 
heavily obscured AGN. We identify a small number of input parameters to the XRB synthesis code that encapsulate 
the minimum level of uncertainty in reconstructing the XRB spectrum. These are the power-law index and high-energy 
cutoff of the intrinsic X-ray spectra of AGN, the level of the reflection component in AGN spectra, and the fraction of 
Compton-thick AGN in the Universe. We then map the volume of the space allowed to these parameters by current 
observational determinations of the XRB spectrum in the range 3-100 keV. One of the least-constrained parameters is 
the fraction of Compton-thick AGN. Statistically acceptable fits to the XRB spectrum at the 68% confidence level can be 
obtained for Compton-thick AGN fractions in the range 5-50%. This is because of degeneracies among input parameters 
to the XRB synthesis code and uncertainties in the modelling of AGN spectra (e.g. level of reflection fraction). The most 
promising route for constraining the fraction of Compton-thick AGN in the Universe is via the direct detection of those 
sources in high-energy (i> lOkeV) surveys. We show that the observed fraction of Compton-thick sources identified in 
the SWIFT/BAT serendipitous survey limits the intrinsic fraction of Compton-thick AGN, at least at low redshift, to 
10-20% (68% confidence level). We also make predictions on the number density of Compton-thick sources that current 
and future X-ray missions are expected to discover. Testing those predictions with data will place tight constraints on 
the intrinsic fraction of Compton-thick AGN as a function of redshift. 

Key words. X-rays: diffuse background; X-rays: general 



1. Introduction 

The origin of the diffuse X-ray background (XRB) has been 
at the centre of high-energy astro physical research since 



its discovery (Giacconi et al. 1962). Already by the early 



k Woltjer 


1989 Awaki et al. 


1991 


). Progress in the inter- 


pretation of the X-ray spectra of AGN (e.g 


Williams et al. 


1992;|Nandra & Pounds 1992; Turner et ai.||1993| |JNandra 


et al.|1993 


Nandra & Pounds 1994 


coupled with improved 



90s numerous X-ray experiments have constrained the most 
salient spectral features of the XRB, the power- law spec- 
tral energy distribution with photon index T = 1.4 in the 
energy interval 2-10 keV and the peak at 20-30 keV (e.g. 
|Marshall et al.|1980|[Gruber et al.|1999[ ). At the same time, 
observations at energies below about 10 keV demonstrated 
that at least a fraction of the XRB intensity is composed 
of point sources, some of which were i dentified with broad - 
line QSOs (Quasi-Stellar Objects, e.g. Shanks et al.||1991 ). 
However, early suggestions that the XRB is the result of 
the superposition of AGN (Active Galactic Nuclei) across 
cosmic time stumbled upon the steeper X-ray spectra of 
UV-bright (unobscu red) QSOs compared to th e XRB spec- 
trum below lOkeV ( |Fabian fc Barcons||T99"2"| ) . This led to 
suggestions that Seyfert-2 AGN, which were observed to 
have flatter X-ray spectra compared to UV bright QSOs, 
could potentially be major contributors to the XRB (Setti 



constraints on their cosmoiogical evolution (Stocke et al 
199Tj paccacaro et al.][l"99l) |Shanks et al.||199H |BcyIe 



et al 



1993) led to the seminal work of Comastri et al 



(19951. They demonstrated that AGN under the unifica- 



tion paradigm (Antonucci 1993) could reproduce the XRB 



spectrum up to about 100 keV. A key prediction of that 
model was the high incidence of obscured accreting super- 
massive black holes (SMBHs) in the Universe, including a 
fraction (w 25% of total population) of luminous Compton- 
thick ones (N H > 10 24 cuT 2 ). 

Observations with the XMM-Newton and Chandra con- 
firmed some of these p redictions, lending supp ort to the ba- 
sic assumptions of the Comastri et al. ( 1995[ ) model. Deep 
surveys carried out by those telescopes resolved almost the 
entire XRB i ntensity into point sources up to energies of 
few keV (e.g. Worsley et al.|2005 ). Moreover, follow-up ob- 
servations suggested that the majority of those sources are 
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AGN (e.g. Bauer et al.|2004 |. Analysis of their X-ray spec- 
tra also demonstrated that in most of them the direct X-ray 
emission is reprocessed by gas an d dust clouds with column 
densities N H « 10 22 - 10 24 cm" 2 ( |Tozzi et al.|2006| | Akylas 



et~aI1[2006l) . X -ray surveys by those telescopes revealed an 



unexpected dependence of the obscured (Nh > 10 cm 2 ) 
AGN fraction on accretion luminosity, wh ich implies devia- 



tions from the basic unification paradigm ( Ueda et al. 
La Franca et al.||2005| |Akylas et al.||2006 



Dwelly et al.||2005| |Dwelly fc Page 



H 



asmger 



2006). 



2003 



2008 



but see 

In the wake of those developm ents, |Gilli et al.| (120071) 
revised the |Comastri et al. ( |1995[ ) XRB synthesis model 
by including the knowledge accumulated by the XMM- 
Newton and Chandra observations both on the evolution 
of the X-ray luminosity function and the dependence of the 
obscured AGN fraction on luminosity. Large numbers of 
Compton-thick sources with column densities 10 24-26 cm~ 2 
(equal number to that of Compton-thin AGN) were re- 
quired by the revised model. Most of them were predicted 
to lie close to, or below, the sensitivity limits of the deep- 
est current X-ray surveys. Attempts to id entify this hidden 
population either directly at X-rays (e.g. Tozzi et al~f2 006 ; 
Georgantopoulos et al.|2009^ |Comastri et al.|2011||Fer uglio 



et al.|20l7 1 or the mid-infrared (e.g. Fiore et al.|2009 ) sug- 
ges ted agreement, at least to the first approximation, with 
the |Gilli et al. (2007) predictions. At the same time, how- 
ever, criticism has been raised on the large percentage of 
Compton-thick AGN adopted by that model. Studies select- 
ing local AGN either at hard X-ray energies <; 15 keV ( e.g 



Beckmann et al.|2009[|Tueller|201l1|Burron et al.|2011|), or 



via o ptical spectroscopy (e.g. |Akylas fc Georgantopoulos 



2009 



2011 



or at the mid-infrared (e.g. Brightman & Nandra 
determined intrinsic Compton-thick fractions in the 
range 15-25% of the AGN population. 

It has also been pointed out that observational uncer- 
tainties as well as aliases among the large number of param- 
eters that are unavoidably adopted in XRB synthesis codes, 
make it hard to draw robust conclusions on the space den- 
sity of Compton-thick AGN in the Universe. |Treister et al. 



(2009), for example, showed that the current level of uncer- 



tainty in the normalisation of the XRB intensity limits the 
predictive power of XRB synthesis codes on the fraction 
of Compton-thick AGN. Moreover, certain AGN spectral 
components are not yet well understood, e.g the strength 
of the radiation reflected off the torus or the accretion disk 
in the line of sight. Tuning this parameter in XRB synthe- 
sis models to values allowed by current observational con- 
straints can change the requi r ed number of Compton-thick 
sources ( |Treister et aT|2009| ) . |Draper fc Ballantynel ( |2009[ ) 
also highlighted the significance of unobscured AGN popu- 
lations for the interpretation of the hard (> lOkeV) X-ray 
background. These authors estimated the contribution to 
the XRB of blazars with X-ray spectra dominated by syn- 
chrotron radiation. Including those sources can reduce the 



fractio n of Compton-thick QSOs required by the Gilli et al. 
(20071) model by about 20%. 



In light of those studies, it is interesting to revisit the 
question of whether the XRB hump at 20-30 keV requires 
the presence of Compton-thick AGN. This paper takes a 
critical approach to this problem by presenting an XRB 
synthesis code that uses state-of-the-art numer ical simula- 
tions to model t he high-energy spectra of AGN ( |Brightman| 
fc Nandra|2011 1. A feature of those simulations is the real- 
istic modelling of the Compton scattering of photons in a 



dense medium, which has a strong impact on the X-ray 
spectra of AGN at obscurations close to and above the 
Compton-thick limit, Nh = 10 24 cm~ 2 . We then identify 
four input parameters to the XRB synthesis code (one of 
which is the number of Compton-thick AGN, others re- 
late to the AGN spectral characteristics), which encapsu- 
late the minimum level of uncertainty in reconstructing the 
XRB spectrum. Instead of fixing those parameters to lit- 
erature values, we allow them to vary within plausible in- 
tervals. We then map the volume of the space allowed to 
those parameters by current observational determinations 
of the XRB intensity. This approach allows us to investi- 
gate in a quantitative way degeneracies among the AGN 
spectral parameters and their impact on the determination 
of the Compton-thick AGN fraction from measurements 
of the XRB spectrum. Finally, we demonstrate how cur- 
rent data and future X-ray missions can provide significant 
progress in our understanding of the origin and composition 
of XRB, by directly constraining the fraction of Compton- 
thick AGN across cosmic time. In our calculations we adopt 
H = 70kms- 1 Mpc- 1 , VL M = 0.3 and Q A = 0.7. 



2. The X-ray background model 

The synthesis of the XRB spectrum requires assumptions to 
be made on (i) the shape of the spectral energy distribution 
of AGN at X-rays, (ii) the volume density of AGN in the 
Universe as a function of accretion luminosity and redshift, 
(iii) the distribution of AGN in obscuration, which is often 
parametrised by the hydrogen column density, Njj- Each 
component of the XRB synthesis model is described below. 



2.1. AGN X-ray spectral energy distribution 

Active Galactic Nuclei have complex spectra at high en- 
ergies that depend on the properties of the SMBH (mass, 
spin), the rate at which matter accretes onto it, and on the 
overall structure and distribution of material in its vicinity 
(e.g. torus). Observations and theoretical work suggest that 
the components needed to reconstruct the X-ray spectra of 
AGN include (i) a power-law with an exponential cutoff at 
high energies, which represents the intrinsic AGN radiation, 
(ii) photoelectric absorption and Compton scattering, (iii) 
emission lines, the most prominent of which is the FeKa 
at 6.4 keV, (iv) reflection of the direct AGN radiation on 
an optical thick medium, (v) an excess emission over the 
power-law extrapolation at soft energies (< 1 — 2 keV), and 
(vi) a possible contribution from radio jets. 

The intrinsic AGN X-ray spectrum can be parame- 
terised by a power-law with photon index T and an ex- 
ponential decline at high energies with e-folding parame- 
ter Eq. This functional form is the result of comptonisa- 
tion processes. Low-energy photons (UV, soft X-rays), pos- 
sibly associated with the accretion disk, are upscattered 
to hard X-ray energies by the energetic particles of a hot 
coro na that surrounds the inn ermost regions of SMBHs 
(e.g. 



Haardt & Maraschi 1991). In this standard model, 
the photon index 1' is related to the temperature, the op- 
tical depth and ove rall structure of the matte r close to the 
central engine (e.g. Murphy & Yaqoob 2011). The energy 
cut-off depends on the temperature of the thermal hot elec- 
trons that are responsible for the Compton-scattering of the 
accretion disk seed photons. 



2 



A. Akylas et al.: X-ray background models 



The direct X-ray emission of AGN is modified by cold 
material in the vicinity of the SMBH, such as the puta- 
tive torus of gas and dust clouds. X-ray photons that in- 
tersect such clouds are absorbed with an efficiency that 
depends on the density and the chemical composition of 
obscuring screen. When the column density of the absorb- 
ing screen increases, Compton-scattcring of X-ray photons 
on the electrons of the intervening material becomes in- 
creasingly important relative to photoelectric absorption. 
The net result is that the direct emission, at least below 
about 10 keV, is significantly suppressed. The Compton- 
scattering optical depth is defined r ~ ijtNh , where ctt is 
the Thomson-scattering cross section and Nh the hydro- 
gen column density of the material. Compton-scattering 
dominates the fate of the X-ray photons emitted by the 
central engine at r ^ 1 or equivalently Nh = (Jt = 
1.5 x 10 24 cm~ 2 . The impact of photoelectric absorption and 
Compton-scattering on the intrinsic AGN spectrum is ac- 
counted for i n a self-consistent way using the numerical 
simulations of Brightman & Nandra (2011). We used their 
model that describes the transmission spectrum from an 
isotropic source located at the centre of a uniform spheri- 
cal distribution of matter. Input parameters of the simula- 
tions are the T and Ec of the intrinsic AGN spectrum, the 
column density Nh, and metal abundance of the obscur- 
ing material. For simplicity we fixed abundances to solar 
values. 

Observations also show the presence of a prominent 
iron Ka line in the spectra of AGN at an energy of 
6.4 keV (iNandra & Pounds! 119941). This can be produced 
as direct AGN emission intersects optically thick material, 
such as the accretion disk or the torus, and is Compton- 
backscattered into the line of sight. The same processes pro- 



relatively soft energies (e.g. 0.5-2 keV). For this application 
we included a soft-excess scattered-light component in the 
spectra of type-2 AGN. It follows the same functional form 
as the direct AGN radiation (power-law plus high-energy 
cutoff) with a normalisation fixed to 3 per cent of the in- 
trinsic AGN spectrum in the 2-10 keV band. This is the 
same as the strength of the soft-excess co mponent for the 
type-2 AGN adopted by |Gilli et ah] p007l ). 

About 10% of the AGN have powerful radio jets that 
may also contribute to or even dominate the emission at 
X-ray wavelengths. Although we did not directly include 
such a spec t ral co mponent, we used the results of |Draper fc 
Ballantyne ( 2009 ) to account for the contribution of blazars 
to the XRB spectrum. 

Given our assumptions for the main emission compo- 
nents of AGN, the key parameters for modelling their in- 
trinsic X-ray spectra are T, Ec, and the strength of the re- 
flection component, Jref- X-ray spectroscopy of individual 
AGN provides constraints on those parameters. There are, 
however, degeneracies among the different emission com- 
ponents which when combined with observational limita- 
tions, do not allow a description of the typical spectral 
characteristics of AGN at the accuracy level required for 
the synthesis of the XRB. For example, the average photon 
index of AGN is determined to be in the range 1.8-2.0 with 
a standard deviation around the mean o f about 0.15-0.30 
( Nandra & Pounds||1994| |Dadina||2008| [Beckmann et al. 
20091 I Brightman fc Jvjandra||2011||Ricci et al.||2011| |Burlon 



et al. 12011 ). Observational constraints on the cutoff energy 



duces a reflection continuum c omponent ( George & Fabian 
1991 Murphy fc Yaqoob|2011 ), which becomes increasingly 



i mportant at e nergies > 8 keV . The numerical simulations 
of Brightman & Nandra (2011 ) include the Ka emission of 
iron as well as the transitions from several other elements. 
There is no reflection component in the s pherical symmet- 
ric model of Brightman & Nandra (2011 1. We chose to de- 
scribe this compon ent using the simulations of |Magdziarz| 2.2. The AGN X-ray luminosity function 



range from as lo w as 50keV (e.g. Molina et al.||2009[ ) to 
about -300 keV ( |Dadina||2008[ ). The strength of the rerlec- 
tion component also varies among different studies, with 
observational con straints ranging from Jref = 0.01 to 0.1 
(|Ricci et al.|2011 ) (corresponding to a reflection parameter 
R — 0.2 — 2) . A possible dependence of the reflection of the 
source characteristics, i.e. type -1 vs type-2 AGN , may ex- 
plain the wide range of values ( Ricci et aT]|2011 ) , although 
such a trend remains to be established. 



& Zdziarski 



( 1995 ) as implemented in the PEXRAV model 
of XSPEC. We parameterise the reflection emission by its 
strength Jref relative to the direct AGN radiation inte- 
grated in the rest-f rame 2-10 keV band. Other XRB syn 



thesis models (e.g. Gilli et al. 2007 Treister et al. 2009) 



measure the amplitude of the reflection component with 
the parameter R, which denotes the solid angle su btended 
by the reflecting material. For comparison R=0.5 ( |Treister| 
et al. [2009" ) corresponds to Jref — 2. 

Ideally, the modelling of AGN X-ray spectra should also 
include a soft-excess component. For type-1 AGN this is 
suggested to be relativi stically blurred photoi onized disc 
reflection emission (e.g. Crummy et al. 20061. In type-2 
AGN, the soft excess is proposed to be scattered radiation 
in the line of sight that is not absorbed by the torus (e.g. 
Turner et al.|1997[ ). When synthesising the XRB spectrum, 
we chose not to include the above components because we 
are primarily interested in the intensity of the XRB spec- 
trum at ^ 2 keV. The contribution of the soft excess is likely 
to be small at those energies. These components become im- 
portant, however, when using the XRB synthesis code to 
make predictions on the expected number of AGN, partic- 
ularly heavily obscured ones, as a function of X-ray flux at 



For the X-ray luminosity function (XLF) of AGNs and its 
evolution w i th red shift we used the parameterisation of 
Ueda et al.| ( |2003h . They estimated the AGN XLF in the 
rest-frame 2-10 keV energy range using a combination of 
hard-band (> 2keV) surveys conducted with the HEAO- 
1, ASCA, and Chandra missions. Here, we adopted the 
luminosity-dependent density evolution of th e luminosity 
functi on, which provides a better fit to the |Ueda et al.| 
(2003) data. In this parameterisation the cutoff redshift 
after which the evolution of AGN stops increases wi th lu- 
minosity. We preferred to use the |Ueda et al. ( 2003 ) XLF 
over more recent determinations (e.g Aird et a L| |2010p be- 
cause the former include a self-consistent estimation of the 
fraction of Compton-thin AGN (N H = 10 22 - 10 24 cm- 2 ). 



2.3. The distribution of AGN in hydrogen column density 

It is assumed that AGN are distributed uniformly in the 
hydrogen column density logarithmic intervals, log Nh = 
20 - 22 and 22 - 24 (cm~ 2 ). The AGN N H distribution, 
J(Nh), is therefore approximated by a step function with 
a fixed ratio, IZ, between absorbed (Nh — 10 22 — 10 24 cm -2 ) 
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Table 1. Range of parameters used in modelling the XRB 



parameter 


min 


max 


step 


(1) 


(2) 


(3) 


(4) 


r 


1.78 


2.02 


0.02 


E C 


125 keV 


300 keV 


5keV 


Jref 


0.01 


0.08 


0.0025 


fcT 


0.05 


0.50 


0.05 



1: parameter of interest, where T and Ec are the power-law 
index and high-energy cutoff of the intrinsic AGN spectrum, 
fnEF is the reflection component strength and fcT is the 
fraction of Compton-thick sources relative to the obscured 
AGN population; 2: minimum value of parameter; 3: maxi- 
mum value of parameter; 4: step size 



and unabsorbed (Nh < 10 22 cm -2 ) sources. The parameter 
1Z is a function of X-ray luminosity as parameterised by 
Akylas et al. (2006). It increases from 1Z ss 0.4 at Lx(2 — 
lOkeV) « 10 4V eTgs- 1 to 11 » 0.8 at L x {2 ~ lOkeV) 



10 42 ergs 1 . For comparison, the corresponding fractions 
in the Gilli et al. (2007) model are K « 0.45 and 0.75 at 
the above luminosities. The adopted parametrisation for 
the luminosity depend ence of the obscure d AGN fraction is 
als o similar to that of l Ueda et al. ( 2003). We chose to use 
the Akylas et al. (2008) results because they are using X- 
ray spectral analysis, not hardness ratios, to determine the 
column density of individual AGN in their sample. It should 
be emphasised that our results and con clusions would not 
change if we used the Ueda et al. ( 2003 ) functional form for 
the Lx dependence of 'It. 



Compton-thick AGN are not included in the Ueda et al. 



(2003) XLF, or indeed any XLF determination based on ob- 



servations in the 0.5-10 keV spectral band. This is because 
for Nh 10 24 cm~ 2 , the bulk of the direct AGN emis- 
sion at rest-frame energies ^ 10 — 20keV is suppressed. 
These heavily obscured sources are therefore too faint in 
the 0.5-10 keV spectral window of Chandra and XMM- 
Newton. Their identification is therefore hard even in the 
deepest current X-ray survey fields (e.g. |Tozzi et "aLj [2006 



Georgantopoulos et al.||2009 ). As a result, the fraction of 



Compton-thick AGN remains uncertain (see Comastri et 
al. 2011 and references therein). We introduced the pa- 
rameter for, which is defined as the fraction of Compton- 
thick AGN (N H > 10 24 cm~ 2 ) relative to mildly obscured 
Compton-thin sources (Nh = 10 22 — 10 24 cm -2 ). We then 
added this fractio n of Compton-thick AGN to the XLF of 
Ueda et al. p003). It is therefore assumed that Compton- 
thick sources follow the same XLF as mildly obscucred sys- 
tems. 



3. Adopted parameters 

Modelling the XRB involves many parameters, some of 
which are reasonably well-constrained (e.g. the XLF of 
Compton thin AGN at least out to z = 2 — 3), while oth- 
ers have larger uncertainties (e.g. AGN X-ray spectrum). 
We identified four parameters that are important for shap- 
ing the XRB spectrum and for which sufficiently accurate 
observational determinations are still not available. These 
are the photon index, T, and the high-energy cutoff, E C1 
of the intrinsic AGN spectrum, the fraction of the reflec- 
tion component in AGN X-ray spectra, /ref, and the frac- 
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Fig. 1. Red curve: an example of an XRB spectrum pro- 
duced by our synthesis model for parameters T — 1.88, 
E c = 230 keV, f REF = 0.055 and f CT = 0.1. The points 
are the observational constraints on the intensity of the 
XRB at different energies from the studies discussed in the 
text. Comparison of the model curve with the data yields a 
X 2 of 70.8 for 86 degrees of freedom (90 data points minus 
four free parameters). 



tion of Compton-thick sources, fcr, among AGN. Instead 
of fixing those parameters to plausible values, we adopted 
the alternative approach of leaving them free to study the 
region of the four-dimensional parameter space, which is 
consistent with current observational data. These include 
measurements of the XRB spectrum itself and the fraction 
of Compton-thick AGN detected by high-energy (<; 10 keV) 
missions. 

We emphasise that the list of free parameters we chose 
is far from exhaustive. For example, a possible r edshift de- 
pendence of the fraction of obscured AGN (e.g. |La Franca] 
eTl^[2005l |Treister fc L7r^y][2006| [Delia Ceca et al.||2008p , 
which is not included in our modelling, would imprint on 
the predicted XRB spectrum. Similarly, the results of the 
synthesis code depend on the spread of the T distribution, 
which we fixed to 0.15 following the BeppoSAX results of 



Dadina (12008). Another limitation of our modelling is that 
we assumed that all AGN have the same values of E c and 
fnEF, when a distribution of values might be more realistic. 
Nevertheless, by investigating the impact of observational 
uncertainties of a few parameters on the reconstruction of 
the XRB spectrum, one can assess whether observations of 
the diffuse X-ray background intensity can place reliable 
constraints on AGN properties. Table [T] shows the adopted 
range of values for each of the four parameters of interest 
as well as the step size used to sample the parameter space. 

For synthesising the XRB spectrum, the assumed XLF 
and Nh function were integrated in redshift, luminosity, 
and hydrogen column density to estimate the expected 
number of AGN at a given bin of redshift, luminosity, 
and Nh- The integration was carried out in the lumi- 
nosity and column density intervals Lx(2 — lOkeV) = 
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Fig. 2. Left: confidence intervals for the model parameters T, Ec, Iref and fcr used to synthesise the XRB spectrum. 
The (blue) thick continuous line corresponds to the 68 per cent confidence level (A% 2 = 4.72 for four parameters). The 
(red) thin dotted line marks the 95th percentile (A% 2 = 9.70 for four paramete rs). Right: Same as right p anel but 
including to our XRB synthesis results the contribution of blazars, as modelled by Draper & Ballantync (2009). 




Fig. 3. Joint constraints on the parameters T, Ec, Jref an d fcr based on bo th the XRB spectrum and the fraction 
of Compton-thick of AGN in the SWIFT/BAT AGN survey ( |Burlon et aLl|2011| ). Left: The (blue) thick continuous line 
marks the region of the parameter space which is consistent at the 68% confidence level with both the observed XRB 
spectrum and the fraction of SWIFT/BAT Compton-Thick AGN. The (red) thin dotted line marks the region which is 
consistent with those o bservations at the 95% confi dence level. Right: Same as left panel including the contribution of 
blazars, as modeled by Draper & Ballantyne (20091, to the XRB synthesis results. 



10 4U - 10 47 ergs-\ N H = 10 2U - lO^cm"- 2 . The red- 
shift range adopted was z=0-7. Limiting the integration 
to z = 4, to explore the impact of a possible exponen- 
tial decline in the AGN XLF at z ^ 3 - 4 (e.g. |Brusa 



et al. 2009 Civano et al. 2011), changes the estimated in- 



tensity of the XRB spectrum by less than 1%. The step in 



z and Lx is logarithmic and variable in size, taking a min- 
imum value of dlogz = 0.1 at z = 1.5 and dlogLx = 0.3 
at \ogLx(2 — lOkeV) = 43.8 ergs" 1 and increasing at 
lower/higher redshift/luminosities. This significantly im- 
proves the computing performance of the code without loss 
of accuracy because the XRB is dominated by AGN at 
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1 - 2 and L x (2 - lOkcV) « 10 43 - lO^ci-gs^ 1 (e.g. 



0.5 



Aird et al.|2010[ ). The N H step size is fixed to A log N H 
cm _z . 

At each z, Lx and Nh bin the observed shape of the 
AGN spectrum a t energies 3-100 keV was dete rmined using 
the simulations of Brightman & Nandra ( 2011 ) to model the 
effect of photoelectric absorption and Compton-scattering 
for a given intrinsic photon index T and high-energy cutoff 
E c . A Gaussian distribution of photon indices was adopted, 
with a fixed spread of a = 0.15 around the mean. Reflection 
was also added to the AGN spectrum using the pexrav 
model of xspec and assuming a fixed fraction /ref of this 
component relative to the direct AGN emission in the 2- 
10 keV band. The resulting spectrum was then normalised 
to the number of AGN within a given z, Lx and Nh bin 
as determined from the integration of the XLF. Adding 
the normalised spectra from all z, Lx and Nh bins yields 
the XRB intensity as a function of energy for a given set 
of parameters T, E c , /ref and fcr- Figure [T] shows an 
example of a model spectrum produced by our code. 

The standard x 2 statistic was estimated for each set of 
parameters (r, E c , /ref, fcr) by comparing the model 
predictions with observations at discrete energies. In this 
exercise we used measurements of t he in tensity of the 
XRB esti mated by |Revnivtsev et al.| ( |2003[ RXTE/PCA, 
3-20keV),|Churazov et al.| (|20071 INTEGRAL, 5-1 00 keV), 



Ajello et al.| |2008b| BW1FT/EAT, 15-200 keV), IMorett 
"1 SWIFT/XRT 



ct al. 



(2010 



Tiirler et al 



mcnts 
on the 
vals, not 



et al. 1995 



p009t SWIFT/XRT 1.5-7 keV) and _ 
INTEGRAL/IBIS, 20-200 keV). These measure- 
are plotted in Figure [T] Observational studies 
XRB spectrum that provide confidence inter- 

(e 



Markevitch 



dis crete data poin ts (e.g. ASCA; 

Kushino et al.| |2002[ Chandra, 



Luca & 



Gcndreau 
HickoTTZ 



2006| XMM -Newton ^Lumb et al" 



Molendi 
used 



2004 
the 



BeppoSAX, 
analysis. The 



Vecchi ct 
TTceV" 



al 



200 2| |De 
1999J 
alisa- 



were not used m ttie analysis. The ikeV norm, 
tions of the XRB spectrum determined from those stud- 
ies are in the range 9.6 — 12keV keV -1 cm~ 2 s _1 sr^ 1 . 
This should be compared with the IkeV normalisation 
of 12.2 keV keV -1 cm~ 2 s _1 sr _1 determined by the 
SWIFT/XRT, which was used as a constraint in our anal- 
ysis. Also, the XRB measur ements from the HEAO-1 mis- 



sion (|Marshall et al 
mode 



1980 ) were not compared with the 
predictions since there is a controversy over the ab- 



solute flux calibration of that mission, which may result in a 
systematic underestimation of the XRB intensity level (e. g 



Ueda et al.|2003[|Worsley et al.|2005 |Treister fe Urry|2005 1 . 
Indeed, almost all measurements of the XRB intensity at 
IkeV lie above the HEAO-1 estimate by factors of up to 
1.5. Stray-light affecting focusing X-ray telescopes, unlike 
the HEAO-1, may also be responsible for the discrepancy. 
Nevertheless, INTEGRAL and SWIFT/BAT, which have 
similar optics to HEAO-1, also found a higher normalisa- 



et al. 2009). 



tion for the XRB intensity (Churazov et al. 2007 Moretti 



The comparison between model and observations is lim- 
ited to the energy range 3-100 keV. Below about 3keV 
sources other than AGN also contribute to the XRB, 
i.e. X-ray emission from ou r Galaxy, groups, and clus- 



ters (Blanchard et al. 1992). At energies above 100 keV 



processes that are not modelled in this paper, such as 
inverse-Compton scattering of photons in relativistic jets o f 
radio- loud AGN, dominate (see Draper & Ballantync 2009). 
These processes are also expected to have a small contribu- 



tion to the XRB below 100 keV ( |Giommi et al |2006| |AjeiIo 
|et al.|[2009| ). We approximated their impact by adding to 
the XRB spectrum estimated by our synthesis code the 



Blazar contribution as modelled by Draper fc Ballantyne| 
(2009), assuming a duty cycle of 13% for the BL Lacs. In 
that model the blazar contribution to the XRB is 12% in 
the 0.5-2 keV band, 7% in the 2-10 keV band, 9% in the 
15-55 keV band and close to 100% in the MeV region. The 
X 2 statistic was re-estimated for each set of our model pa- 
rameters (r, E c , /ref, }ct), including the contribution of 
blazars. 

The XRB synthesis code also produces predictions on 
the number density of AGN samples that fulfill different se- 
lection criteria in luminosity, redshift, and/or obscuration. 
The observed fraction of Compton-thick AGN as a function 
of flux limit and detection energy band provides constraints 
on the parameters of the XRB synthesis model that are in- 
dependent of those obtained from the XRB spectrum. In 
this exercise we used the fraction of Compton-thick AGN 
detecte d in the 15-55 keV en ergy band by the SWIFT/BAT 
survey (Ajello et al. 2008a). The sky coverage sensitivity 
curve corresponding to the 5 a detection threshold of the 
S WIFT/BAT sample was folded in our XRB synthesis code 
to predict the fraction of Compton-thick sources in that sur- 
vey for the grid of model parameters presented in Table [I] 
This was then compared with the observations using bino- 
mial distribution confidence intervals. There are 197 AGN 
in the SWIFT /BAT s ample of which nin e (4.5%) are clas- 
sified Compton-thick ( |Burlon et al.|2011[ ). The 68 and 95% 
binomial confidence intervals for the Compton-thick frac- 
tion are 0.031-0.060 and 0.016-0.074, respectively 



4. Results 



Figure [2] plots the two-dimensional projections of the four- 
dimensional parameter space, which is consistent with the 
XRB observations. All possible pairs of the parameters T, 
E c , /ref-, and fcr are shown in different panels. The con- 
tours correspond to the 68 and 95% confidence intervals. 
In the case of four degrees of freedom these correspond to 
A\ 2 relative to the minimum x 2 value of 4.72 and 9.70, 
respectively. The right set of panels in Figure [2] incl ude the 
contribution of blaza rs to the XRB as modelled by [Draper 



fc Ballantyhe]p0fJ9 |. 

The parameter space shown in Figure [2] is complex. 
Firstly, there are degeneracies among parameters, such as 
r and Ec, with higher values of the photon index also re- 
quiring higher cutoff energies to provide acceptable fits to 
the XRB observations. There are also effects related to the 
resolution of the parameter grid. These are manifested by 
wiggles in the contours e.g. in the plots of fcr vs Ec or 
/ref vs Ec- Increasing the sampling rate of the parame- 
ters (i.e. smaller step size in Table [T]) would smoothen the 
contours at the expense of computation time. 

Despite this limitation, the XRB spectrum does place 
constraints on some of the four parameters plotted in Figure 
[2] For example, there is a relatively narrow region of the 
T-Ec parameter space that is consistent with the observa- 
tions at the 68 or the 95 % confidence interval. In contrast, 
all plots that include the fcr parameter (top row of panels 
in Figure [2]) have confidence contours that cover most of 
the grid space. Therefore the parameter that is least sensi- 
tive to the XRB spectrum (and hence least constrained) is 
the fraction of Compton-thick AGN. For any value of fcr 
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Table 2. Model parameters (T, Ec, Jref, fcr) used to 
predict the fraction and number counts of Compton-thick 
AGN in Figure [4j These parameters yield XRB spectra 
that are consistent with the observations at the 68% con- 
fidence level. They are also consistent with the frac tion of 
Compton-thi ck AGN in the SWIFT/BAT sample of |Ajello 
et al. (2008a) at the 95% confidence level. 



r 


Ec 


Iref 


fcT 


1.9 


245 


0.0650 


0.05 


1.9 


245 


0.0600 


0.10 


1.9 


240 


0.0575 


0.15 


1.9 


235 


0.0550 


0.20 


1.9 


230 


0.0525 


0.25 



there are combinations of the other three parameters that 
yield acceptable fits to the XRB spectrum at the 68 or 95% 
confidence intervals. 

Including blazars into the X-ray background synthesis 
modelling reduces the volume of the parameter space, which 
is consistent with the XRB spectrum at the 68% confidence 
level. This trend is more pronounced for fcr- Blazars have 
hard X-ray spectra, dominated by synchrotron emission, 
and therefore their net impact on the XRB modelling is the 
reduction of the fraction of Compton-thick sources needed 
to reproduce th e shape and overall norma lisation of the 
XRB spectrum (Draper & Ballantync 2009). Nevertheless, 
blazars appear to have a small impact on the 95% confi- 
dence level contours of Figure 2] which remain wide, espe- 
cially for panels that include the fcr parameter. 

Figure [3] shows how the (T, Ec, Iref, fcr) parame- 
ter space is modified by considering the joint constraints 
from both the XRB spectrum and the observed fra ction of 
Compton-th ick AGN in the SWIFT/BAT survey ( |Burlon 
et al.|[20lT ). As expected, the most affected parameter is 
fcr, which is limited to the range 10-20% (68 per cent con- 
fidence level). This is similar to the Compton-thick AGN 
fraction i n the local Universe estimated by optical and IR 
surveys fl Akylas k, Georgantopoulos 2009 Brightman & 
Nandra||2011| ). 



5. Predictions on the properties of Compton-thick 
AGN 



Current and future 
eROSITA 



2010 



e.g. NUST AR, |Harrison et al i 
Predehl et al.||2011|) high ssions have 

among their prime scientific goals the direct identification 
of Compton-thick AGN to constrain their cosmological evo- 
lution. This section presents the predictions of our model on 
the number density and redshift distribution of Compton- 
thick AGN at different energy bands and depths. 

Figure [4] plots as a function of flux the observed fraction 
and cumulative number counts of Compton-thick sources 
determined by our XRB synthesis code for different values 
of the parameter fcr- The curves in that figure were con- 
structed by fixing r — 1.9 and then setting E cut , Jref to 
the values that yield the minimum \ 2 f° r the XRB spec- 
trum (see Table [2]). Predictions are presented in the 0.5-2, 
0.5-8, 2-10 and 20-40 keV bands. The redshift distributions 
of the Compton-thick AGN at different flux limits in those 
bands are shown in Figure [5] The curves in that figure cor- 
respond to the model of Table [2] with fcr = 15%. 



The 0.5-2 keV band is not an obvious choice of en- 
ergy interval for Compton-thick AGN searches. However, 
many current (e.g XMM-Newton, Chandra) and future (e.g. 
eROSITA) missions have their highest throughput in that 
energy range. Therefore, the deepest current images of the 
X-ray sky, in terms of resolving the highest fraction of 
the XRB, are in the 0.5-2 keV band. The dctcctability of 
Compton-thick AGN in the 0.5-2 keV band depends on both 
the assumed fcr and the adopted level of the soft compo- 
nent (e.g. scattered radiation) of the AGN spectra. The lat- 
ter parameter has a strong impact on the number density 
of Compton-thick AGN predicted by our model. 

The choice of the 0.5-8 keV band is motivated by recent 
results on the numb er density of Compton-th ick AGN in 
that spectral band by Brightman & Ueda ( 2012 ). They used 
the 4 Ms Chandra Deep Fie ld South data in comb i natio n 
with the spectral models of Brightman fc Nandra| ( |2011[ ), 
which account for Compton-scattering and the geometry 
of the circumnuclear material, to identify Compton-thick 
AGN to the flux limit f x (0.5-8 keV) « 10" 16 ergs^ 1 cm" 2 . 
Figure [4] shows that the cumulative number counts of those 
sources are consistent with the fcr = 25% model curve. 
We chose to compare our model predictions with the ob- 
served cumulative number count distribution of Compton- 
thick sources rather than their fraction relative to the over- 
all AGN population. This is to avoid any uncertainties in 
the determination of the number density as a function of 
flux of the overall AGN population with Lx ( 2 — l OkeV) > 
10 42 ergs _1 adopted by Brightman & Ueda (2012) to nor- 



malise the Compton-thick number counts. In the context 
of our model the f ractio n fcr ~ 25% implied by the 

h 



Brightman & Ueda (|2012 ) results is consistent at the 95% 
confidence level with the observed fraction of Compton- 
thick sources in the SWIFT/BAT AGN sample in the local 
Universe (see Figure [3]). This tentatively suggests a redshift 
dependence of the fcr, given the mean ex pected redshift, 
z « 1, of the |Brightman_fe Ueda (2012) Compton-thick 
AGN sample (see Figure [5]). 

The 2-10 keV band is routinely used to search for heav- 
ily obscured AGN. This is because photons at this energy 
interval are least affected by moderate obscuring column 
densities, while the sensitivity of imaging telescopes, such 
as XMM-Newton and Chandra, at this energy range re- 
mains relatively high. The detectability of Compton-thick 
AGN in the 2-10 keV band depends on the level of the re- 
flection component of the AGN spectra. 

The 20-40 keV band encompasses the peak of the XRB 
spectrum and will be explored to unprecedented depths 



10- 



ergs cm ) with the imaging optics of the 



NUSTAR mission (Nu clear Spectroscopic Telescope Array, 
Harrison et al.||2010"| . The Compton-thick AGN fraction 
predictions in that band are more directly related to fcr 



predictions in that band are more directly related to fcr 
and are least dependent on other input parameters to 
the XRB synthesis code, such as the level of the reflec- 
tion component. At the confusion flux limit of NUSTAR, 
/x(20 — 40keV) « 10 -14 ergs -1 cm -2 , we predict a frac- 
tion of Compton-thick AGN between 2 to 10%, depending 
on fcr- Figure [5] shows that the redshift distribution of 
those sources is expected to peak at z ~ 0.6. 

Although the efficiency of finding Compton-thick AGN 
is highest in the 20-40 keV band, the largest numbers of 
such sources will be produced by wide-area surveys at lower 
energies. At the flux limit /x(0.5 — 2) « 10 _14 ergs _1 cm -2 
of the eROSITA All Sky Survey for example, the predicted 
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Fig. 4. The observed fraction of Compton-thick AGN in the spectral bands 0.5-2, 0.5-8, 2-10 and 20-40 keV. Our model 
predictions for an intrinsic Compton-thick AGN fraction of 5, 15, 25% correspond to the the red dotted, blue short-dashed 
and gree n long dashed cur ves, respectively. The black solid lines are the prediction s of the X-ray backgrou nd synthesis 



model of Gilli et al. (20071. The magenta dot-dashed lines are the predictions of the Ballantyne et al. ( 201lj ) model. The 
inset plots show the predictions for the cumulative number counts of Compton-thick AGN . The curves are the same as in 



the m ain panels. The observ ational datapoi nts are from T ozzi e t al. (120061) (2-10 keV), Brightman & Ueda (20121 (0.5 

Mr 



keV), |Krivonos et al.| p007| ) adapted from [Ireister et al.| ( |2009p (20-40keV) 



number density of Compton-thick AGN is 0.2-0.9deg 2 . 
This translates into a total of 900 to 4000 Compton-thick 
AGN in the 0.5-2 keV eROSITA source catalogue. Figure[5] 
shows that the mean expected redshift of these sources is 
z « 0.3. In the 2-10 keV band, eROSITA will achieve a sen- 
sitivity of about 3 x 10 _13 crgs _1 cm~ 2 . A total of about 
70 to 280 Compton-thick AGN are expected at that flux 
limit over the whole sky. The relatively small field of view of 
NJJSTAR (6 arcmin radius) and the fact that it is confusion 
limited at / x (20- 40) « 10~ 14 erg s" 1 cm' 2 (about 100 ks 
exposure), means that this mission will detect up to a few 



hundred Compton-thick AGN, depending on the adopted 
survey strategy. A limitation of the eROSITA data, how- 
ever, is that the detected sources will typically have limited 
number of counts. The identification of the Compton-thick 
AGN among the eROSITA X-ray source population will 
therefore be challenging. 



6. Discussion 

We presented a synthesis model for the XRB that uses the 
results of state-of-the-art Monte Carlo simulations for the 



8 



A. Akylas et al.: X-ray background models 



1 1 I I -H F 1 1 I I I I I | 1 1 I I I I llj 1 1 I I I I II 1 1 I H 




0.001 0.01 0.1 1 0.001 0.01 0.1 1 

redshift redshift 



T 1 I I I I 1 1| 1 1 I I I 1 1 1 1 1 1 I I I I 1 1| 1 1 I I -H F 1 | 1 | 1 | 1 1 I H 




redshift redshift 

Fig. 5. The differential redshift distribution of Compton-thick AGN (dN/dz) predicted by our model in the spectral 
bands 0.5-2, 0.5-8, 2-10 and 20-40 keV. The different curves correspond to flux limits of 10~ 13 (black dash-dotted), 10~ 14 
(blue large-dashed), 10 -15 (cyan small-dashed), 10 -16 (red dotted) and 10 -17 ergs -1 cm -2 (red continuous). 



key spectral components of AGN. This was used to recon- 
struct the XRB spectrum for a wide range of values of the 
parameters (i) T, the exponent of the intrinsic AGN power- 
law, (ii) Ec, the high-energy cutoff of the AGN spectra, 
(iii) fiiEF, the level of reflection in the spectra of AGN, 
and (iv) fcr, the fraction of Compton thick sources. We 
then define the volume of the four-dimensional parameter 
space for which the model XRB spectra are consistent with 
the observed one. 

A result from Figures [2] and [3] is that there are no 
solutions consistent with the observed XRB spectrum for 
T 1.98 (better than 95 per cent confidence level), under 
the assumptions of the model and for the range of param- 
eter values explored. Therefore, within the context of the 
model, this is an upper limit to the intrinsic power-law in- 
dex of AGN. 



It is also interesting in Figure [2] that the intensity of the 
integrated XRB spectrum does not provide any meaningful 
constraints on the space density of Compton-thick AGN in 
the Universe. Our analysis shows that a wide range of val- 
ues of the Compton-thick fraction, fcr, yields acceptable 
fits to the XRB spectrum. This is because of degeneracies 
among input parameters to the XRB synthesis code. Also 
the parameters that describe the shape of the AGN spec- 
tra (e.g. reflection fraction, high energy cutoff) are still not 
adequately constrained by observations and therefore the 
range of plausible values remains wide. Figure [3] shows that 
the most solid limits on the fraction of Compton-thick AGN 
comes from the direct detection of such sources. 

This is also demonstrated in Figure |4j which plots 
the model expectations for the observed fraction and the 
number-count distribution of Compton thick AGN as a 
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function of flux at different energy bands. The number 
counts of Compton thick AGN in the Chandra Deep Field 
South (Brightman & Ueda 2012) suggest fcr ~ 25 per cent 
at z w 1, in agreement at the 95 per cent confidence level 
with the Swift/BAT results in the local Universe. 

Figure [4] also compares our predictions on the num- 
ber density of Compton-th ick AGN w i th th e X-ra y back- 



ground synthesis m odels of |Gilli et al. ([2007 Fj and [Draper 

"lit 



& Ballantyne ( 2010 ). The former model is built on the same 
basic principles (e.g. AGN unification paradigm) as the 
XRB code presented in this paper. There are important dif- 
ferences between the two codes however, in e.g. the adopted 
X-ray luminosity function, the overall fraction of Compton- 
thick AGN, and the level of the reflection component. Only 
mildly Comp ton-thick AGN, l ogN H = 24-25 cm -2 , are in- 
cluded in the Gilli et al. ( 2007 ) predictions plotted in Figure 
[4] to facilitate the comparison with our results. Cont r ary to 
other XRB synthesis codes, Draper & Ballantyne ( 2010 ) 
postulated that the physical parameter that is more rele- 
vant to obscuration and Compton-thickness is the accre- 
tion rate onto the SMBH. This is motivated by numerical 
simulations, which suggest that the AGN obscuration is re- 
lated to the evolutionary status of the SMBH, and by obser- 
vational constraints on the accretion properties of nearby 
Compton-thick QSOs. This model assumes that the frac- 
tion of Compton-thick AGN (log Njj = 24 - 25 cm" 2 ) is a 
strong function of accretion rate. They dominate at both 
high (A > 0.9; for « 0.86) and low (A < 0.01; fcr » 0.60) 
Eddington ratios, while they are absent at intermediate 
Eddington ratios. 

In all panels of Figure |4] the IGilli et al.l p007| and 
Draper & Ballantyne (20101 results predict a very diffcr- 
ent flux dependence of the observed Compton-thick frac- 
tion and number count distribution. In the 0.5-8 or 2- 
10 keV bands for example, the differences are caused by 
the adopted fcr and the treatment of the reflection spec- 
tral component in the different codes. We a s signed the same 



Jref to all AGN, whereas in Gilli et al. (20071 the level 
of reflection depends on the obscuration a nd accretion lu- 
minosity. As a result, the IGilli et al.|(|2007|) predictions cut 



through our curves in Figure [4[ At faint fluxes, below about 
10" 16 ergs" 1 cm" 2 , their curve lies well above our predic- 
tions because of the higher fra ction of Compton-thick AGN 
adopted by their model. The Draper & Ballantyne (20101 
model predicts more Compton-thick AGN relative to the 
Gilli et al. (2007) and our predictions. The hump in their 
number-count predictions marks the change in the domi- 
nant Compton-thick population from high Eddington ratio 
sources at bright fluxes to slowly accreting systems. 

The detectability of Compton-thick AGN in the 0.5- 
2 keV band depends quite strongly on the level of the soft 
spectral component. It is not surprising that in Figure [4] 
there is fai r agreement betw een our fcr — 25% model 
and that of [Gilli et aL] p007[ ) down to f x {0-5 - 2keV) « 
10" 16 ergs cm -z , although the latter assumes a much 
higher fraction of Compton-thick AGN. This is because 
both models use the same level of scattered radiation (3%) 
to model the soft component of AGN. The upturn in the 
fraction of the Compton-thick AGN at /x(0.5 - 2keV) ^ 



10" ib ergs" 



in the Gilli et al. (2007) models is caused 



by the increasing contribution of high rcdshift sources for 



which the reflection component is rcdshiftcd into the 0.5- 
2 keV band. 

Finally, the detection rate of Compton-thick AGN in the 
20-40 keV band is more directly related to the total number 
of sources in the parent population and less sensitive to the 
spectral characteristics of AGN (e.g. level of scattered radi- 
ation or reflection fraction) . Our estimates for the observed 
fraction of Compton-thick AGN in that ban d ar e therefore 
consistently below those of | Gilli et al. ( 2007 1 and Draper & 



Ballantyne (2010) because these aut 



http: / /www. bo. astro.it/~gilli/counts. html 



lors assume a higher 
intrinsic fraction of Compton-thick AGN. 

7. Conclusions 

A model for the synthesis of the XRB was presented, which 
used as input an AGN X-ray spectral library generated 
by Monte Carlo simulations to account for the impact of 
Compton-scattering and photoelectric absorption. We then 
explored whether the input parameters to the XRB syn- 
thesis code can be constrained by observations of the XRB 
intensity as a function of energy. Four parameters were used 
in this exercise, the power-law index and high-energy cut- 
off of the intrinsic X-ray spectrum of AGN, the level of 
the reflection component in AGN spectra, and the intrinsic 
fraction of Compton-thick AGN. We mapped the volume 
of this four-dimensional space that is consistent with ob- 
servational determinations of the XRB spectrum in the en- 
ergy interval 3-100 keV and found that one of the least con- 
strained model parameters is the fraction of Compton-thick 
AGN in the Universe. Acceptable fits to the XRB spectrum 
(68% confidence level) can be obtained for Compton-thick 
AGN fractions in the range 5-50%. This is because of aliases 
among input parameters to the XRB synthesis code as well 
as because of uncertainties in our current understanding 
of AGN spectra (e.g. level of reflection). The most strin- 
gent constraints on the fraction of Compton-thick AGN 
in the Universe comes from the direct detection of such 
sources in high-energy surveys, such as those carried out by 
SWIFT/BAT and INTEGRAL. These observations suggest 
a Compton-thick fraction, at least in the local Universe, of 
10-20% at the 68% confidence level. We also predicted the 
number counts of Compton-thick AGN as a function of flux 
and energy band. These predictions can be tested by cur- 
rent and upcoming X-ray missions to set tight constraints 
on the space density of Compton-thick AGN as a function 
of redshift. 
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